arXiv:1507.03357v3 [physics.atom-ph] 11 Oct 2016 


Quantum mechanical calculation of Rydberg-Rydberg autoionization rates 


Martin KifTner^’^, Davide Ceresoli^, Wenhui and Dieter Jaksch^’^ 

Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 

Clarendon Laboratory, University of Oxford, Parks Road, Oxford 0X1 3PU, United Kingdom^ 
Istituto di Scienze e Tecnologie Molecolari CNR, via Colgi 19, 20133 Milano, Raly^ and 
Department of Physics, National University of Singapore, 117542, Singapore"^ 

We present quantum mechanical calculations of autoionization rates for two Rubidium Rydberg 
atoms with weakly overlapping electron clouds. We neglect exchange effects and consider tensor 
products of independent atom states forming an approximate basis of the two-electron state space. 
We consider large sets of two-atom states with randomly chosen quantum numbers and hnd that 
the charge overlap between the two Rydberg electrons allows one to characterise the magnitude of 
the autoionization rates. If the electron clouds overlap by more than one percent, the autoionization 
rates increase approximately exponentially with the charge overlap. This Ending is independent of 
the energy of the initial state. 
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I. INTRODUCTION 

Exciting ultracold atoms to Rydberg states [T] with 
large principal quantum number n furnishes the atoms 
with extremely exaggerated properties. For example, 
the size, interaction strength and polarizability increases 
by several orders of magnitude as compared to ground 
state atoms. This feature allows one to study funda¬ 
mental physical phenomena on completely new time and 
length scales and magnifies physical effects such that they 
become experimentally accessible. For example, dipole- 
dipole interactions between ground state atoms are typi¬ 
cally weak, but they are strong and long-ranged between 
Rydberg atoms such that ^m-sized molecules consisting 
of two w and three [TOHH] atoms become possible. 
Moreover, dipole-dipole interactions between Rydberg 
atoms give rise to the blockade effect [laiii] and crys¬ 
tals of spatially ordered Rydberg excitations that were 
experimentally observed in m- The modification of 
the quantum dynamics of Rydberg electrons due to their 
dipole-dipole interaction has been demonstrated in a re¬ 
cent experiment by Takei et al. m via ultrafast pump- 
probe laser techniques. 

In systems of dipole-dipole interacting Rydberg atoms 
the interatomic spacing is typically large compared to 
the size of the Rydberg electron orbital. A fascinat¬ 
ing prospect for future studies is the investigation of 
Rydberg systems with overlapping electron clouds. In 
this regime the exaggerated properties of Rydberg atoms 
would allow one to study the rich physics of electron- 
electron interactions on much more accessible time and 
length scales compared to conventional solid state sys¬ 
tems. More specifically, the size of the valence electron 
orbital increases like i?„ = 4n^ao where oq is the Bohr 
radius and will thus reach the typical separation between 
atoms in optical lattices or tweezers for n ^ 35. Overlap¬ 
ping electron clouds could give rise to delocalized elec¬ 
trons and correlated quantum many-body states via the 
strong Coulomb interaction between the electrons. How¬ 
ever, these coherent processes compete with autoioniza¬ 


tion and radiative decay processes enabled by the large 
number of empty orbitals below the Rydberg state. 

The first step in investigating the regime of Rydberg 
atoms with overlapping electron clouds is to characterize 
the time scales of the occurring physical processes. While 
radiative processes are well understood m, we here fo¬ 
cus on autoionization of two neutral Rydberg atoms via 
the Penning effect [T] as shown in Fig. In this process 
the energy for ionizing atom A is provided by a change 
in internal energy of atom B. Until now, quantum me¬ 
chanical calculations of this process are restricted to the 
dipole-dipole interaction regime of non-overlapping elec¬ 
tron clouds HHms] where the decay rates are negligibly 
small. On the contrary, calculations based on classical 
Hamilton equations [20] show that fast autoionization oc¬ 
curs for atomic separations of the order of i?„ where the 
electron clouds start to overlap. This effect has been 
identified in [5T] as a key factor for understanding the 
fast autoionization of a Rydberg gas observed in [22] . In 
order to determine the timescale of autoionization of Ry¬ 
dberg atoms with overlapping electron clouds, quantum 
mechanical calculations of the corresponding autoioniza¬ 
tion rates are needed. However, a rigorous approach to 
this problem is extremely challenging since it would in¬ 
volve finding the highly excited two-electron eigenstates 
of the system. 

In order to estimate autoionization rates of Rydberg 
atoms with overlapping electron clouds, we here present a 
simplified model and consider two-atom states IV'm) with 
weakly overlapping electron clouds as shown in Fig. 
We assume that \iPm) is a tensor product of two generally 
different independent-atom orbitals and neglect exchange 
effects. In order to account for the fact that these states 
are not eigenstates of the system, we consider large sets of 
states IV’m) with different quantum numbers that could 
serve as an approximate basis of the true two-electron 
eigenstate. We evaluate the autoionization rate of the 
states IV’m) quantum mechanically and show that the 
charge overlap between the two atoms allows one to char¬ 
acterize the magnitude of the autoionization rates. In the 
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regime of very small charge overlap between the Rydberg 
orbitals, the autoionization rates are small and depend on 
the energy of the initial state. Moreover, the full interac¬ 
tion Hamiltonian can be approximated by its multipole 
expansion. On the contrary, above a certain threshold 
the multipole expansion becomes invalid and the autoion¬ 
ization rates increase approximately exponentially with 
the charge overlap. 

Note that the autoionization mechanism between two 
Rydberg atoms considered here is related to autoioniza¬ 
tion processes in crystals and clusters that have been 
termed inter-atomic Auger decay |23j and more recently 
interatomic Coulombic decay (ICD) [21142 7j . In par¬ 
ticular, the strong enhancement of autoionization rates 
through the overlap between electron orbitals in clusters 
was reported in [^, and ICD processes between several 
excited atoms in a cluster were studied in [23127] . 

This paper is organized as follows. The system of inter¬ 
est and our model are described in Sec.jTIj We briefly out¬ 
line the calculation of the autoionization rates in Sec. IIIII 
and defer more technical details to Appendix]^ In order 
to account for many different initial states "we ran¬ 
domly select these states as described in Sec. |TVj Finally, 
the autoionization rates of the randomly selected states 
are presented in Sec. and a conclusion of our work is 
given in Sec. lYg 

II. THE SYSTEM 

We consider two Rydberg atoms as shown in Fig. [^a), 
where each atom is comprised of a singly-charged core 
and one valence electron. We assume that the atoms 
are so cold that their positions do not change during the 
decay process. Atom A is centered at the origin and atom 
B is located at i? = i?e^, where is the unit vector in 
z direction and R is the atomic separation. The total 
Hamiltonian of the two-atom system is H = Hq + V , 
where Hq = Hq^^ + Hq^^ and Hq^'' is the Hamiltonian 
of Rydberg atom X. All interactions between atom A 
and atom B are described by 

+ ^_^_(1) 

47r£o \ri-f 2 \ \r 2 \ \R-ri\)' 

where q is the elementary charge and the operator 

associated with the position of electron i. The first term 
in Eq. Q accounts for the repulsion of the two ion cores, 
the second is the electron-electron interaction and the 
third (fourth) term describes the interaction of electron 
2 (1) with ion core A (B). The eigenstates of Hq are 

I'^m) = ( 2 ) 

where and |^b) are independent-atom Rydberg 

wavefunctions centered at the origin and i?, respectively, 

tpA{ri) = l/juAUmAiri), tl}B{r2) = 'lpnBlBmB{r2 - R)- 

( 3 ) 
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FIG. 1. (Color online) (a) Two Rydberg atoms in differ¬ 
ent electronic states and with spatial separation R. The blue 
spheres indicate the size of the Rydberg electron charge den¬ 
sity cloud, and Vi are the electron coordinates, (b) Schematic 
illustration of the autoionization process. Atom 2 makes a 
transition into a lower bound state and the electron of atom 
1 is ejected into the continuum, (c) Charge density p for wave 
function ^'m in Eq. 0 with ua = 20, IA ~ 2, ruA = — 1 and 
ub = 15, Ib = 1, rnB = 1. The dashed sphere with radius Ra 
(Rb) denotes the classical outer turning point of the electron 
of atom A (B). 

We ignore the fine structure such that the wavefunctions 
4’nimiT) are characterised by the principal quantum num¬ 
ber n, the orbital angular momentum quantum number 
I and the azimuthal quantum number m. We generate 
the functions ipnim with energy Eni = —l/[n — 5n{l)Y 
via the Numerov method, where Sni}) is the quantum 
defect [DEH]. We choose Rubidium 85 atoms which are 
a popular choice in recent Rydberg experiments and ob¬ 
tain the energies Eni (and hence the quantum defects) 
for n < 11 from spectroscopic data reported in [29]. For 
n = 11, the non-zero quantum defects are Jn(0) =3.134, 
(5ii(l) = 2.652 and Sii{2) = 1.341 which is consistent 
with the quantum defects provided in m- We ignore 
the weak dependence of Sn{l) on the principal quantum 
number for n > 11. Note that we order the quantum 
numbers in |'0 m) such that E^aIa ^ ^ubIb by conven¬ 
tion since the state obtained by interchanging A and B 
has the same autoionization rate. 

An example for I'i/'M) is shown in Fig. [^c), where the 
size of the electron cloud of atom A (B) is indicated by a 
sphere of radius Ra {Rb), where Ra (Rb) is the classical 
outer turning point, 

Rx = nx + nx \J n*x^ — Ixilx + 1)) 


( 4 ) 
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n* = {n — (5„(?)] is the effective quantum number and 
X G {A,B}. In order to quantify the overlap between 
the wavefunctions \iPa) and |'0 b), we consider the amount 
of charge due to \ipA) inside the sphere Vb with radius 
Rb around atom B, 

SqAiVB)=q J I'lpAl'^d^r. ( 5 ) 

Vb 

Similarly, 

6qBiVA)=q J I'tpBfd^r ( 6 ) 

Va 

is the amount of charge due to |'0 b) inside the sphere 
Va with radius Ra around atom A. A measure for the 
differential overlap between \iPa) and |'0 b) is then given 
by 6qjQ, where Q = 2q is the total electron charge and 

Sq = SqAiVB) + SqBiVA). ( 7 ) 

In the following we will consider only states with small 
overlap such that 5q/Q 1. 

Note that the physical wavefunction of the two-electron 
system should include spin degrees of freedom and be 
completely antisymmetric with respect to electron ex¬ 
change. However, we find that the simplified state in 
Eq. § results in a good approximation of the autoion¬ 
ization rate for weakly overlapping electron clouds as ex¬ 
plained in Sec. |III| 

III. AUTOIONIZATION RATES 

Next we outline the calculation of the autoionization 
rate for state \tpM) in Eq. To this end, we consider 
the process shown in Fig. iRb) where atom B makes a 
transition to a lower bound state ipbi'i' 2 ) = ’ipnbhmti'i — R) 
with energy Eb = and the other electron is 

ejected into the continuum. We model the wavefunc¬ 
tion of the ejected electron with mass me by energy- 
normalized Coulomb waves \tp^) [SH 132] with angular 
momentum Z, magnetic quantum number m and energy 
E obeying the generalized normalization relation 

The Coulomb waves are numerically generated by fol¬ 
lowing the procedure described in [Tj. We calculate the 
autoionization rate using Fermi’s golden rule [33] and to 
first order in the interaction V. The decay rate T^ for 
the process shown in Fig. [^b) is thus given by 

(y ^ k 

= E (9) 

lk=0 mk = -lk 

where the energy = Em — Eb of the Coulomb wave is 
fixed by energy conservation between the initial and final 
states and 

Em = {i’M\H\'ipM) (10) 
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FIG. 2. (Color online) Distribution of the principal quantum 
numbers ua and ns in the randomly chosen sets of states Si. 
(a) Red dots correspond to set <Si. (b) Green squares show 
set 52. (c) Blue triangles are for set S3, (d) Black diamonds 
correspond to set 54 . 


is the expectation value of the total Hamiltonian H 
in the initial state \iPm)- Note that Em differs by at 
most 2% from the unperturbed value ('0 m|77o|^m) = 
EuaU + EubIb for S'!! states considered in Sec. [TV] This 
is consistent with the wavefunctions comprising the ini¬ 
tial state being only weakly perturbed by the electron- 
electron interaction in the overlap region. 

A more rigorous calculation with a fully antisymmet¬ 
ric initial state would result in two Coulomb matrix el¬ 
ements in Eq. ® that are termed the direct and the 
exchange term |33) . The single matrix element in Eq. 
corresponds to the direct term, and the exchange term is 
absent since our initial state |'0 m) in Eq. © is a simple 
product state. The exchange term depends on the over¬ 
lap between the single-electron orbitals and decreases ex¬ 
ponentially with increasing distance R [25j . Since we are 
considering only weakly overlapping electron clouds we 
expect that exchange effects are small and hence the ex¬ 
pression in Eq. ([^ should be a good approximation for 
the autoionization rate. 

The decay rate T^^ in Eq. Q accounts for all processes 
where atom B makes a transition into the bound state 
\ipb) and atom A is ionized. In addition, we consider also 
the autoionization process with rate T^ where atom A 
makes a transition to |'(/'b) and atom B is ionized. The 
full decay rate is then obtained by adding T^ and T^ 
and summing over all bound states, 

rM= ^ (r^ + f^M). (11) 

b with 

Eb<EM 

We numerically evaluate Eq. 0 by restricting the sum 
over bound states to those with nb > 0.2(n^ -bns). This 
is justified since the contribution of lower-lying bound 
states is negligible. The evaluation of the matrix ele¬ 
ment in Eq. (jllj) is described in detail in Appendix 
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In short, we expand all involved wavefunctions and the 
interaction Hamiltonian V in Eq. Q in terms of spher¬ 
ical harmonics and limit the integration region to the 
volume where \iPa,b) both take on non-negligible values. 
We restrict the maximum angular momentum in the ex¬ 
pansion of the wavefunctions to I = 1000, and all terms 
in the expansion of V leading to an exchange of angular 
momentum > 15 between the electrons due to the 
Coulomb interaction are neglected. With these choices 
the numerical expense of calculating one value of Em still 
takes up to 20 hours on a 16 core Intel E5-2640v3 com¬ 
pute node. We estimate that the numerical uncertainty 
in Em due to these approximations is approximately 10% 
for initial states with Sq/Q > 10“^, while we achieve full 
convergence for states with 6 q/Q < 10“^. 


IV. SELECTION OF RANDOM STATES 


The two-atom states in Eq. ([^ are independent-atom 
states and thus not eigenstates of the total Hamiltonian 
H. For sufficiently small atomic separations R, the inter¬ 
action V couples many states |V’m) with different quan¬ 
tum numbers mm- Here we are not interested in the 
quantum dynamics of a particular initial state, but the 
aim is to characterize the autoionization rates of a large 
variety of different states To this end, we calculate 

the autoionization rate Em for four sets Si of randomly 
chosen states that we select as follows. We consider 
four non-overlapping energy intervals that are centered 
around the energies Si, S 2 , S 3 and S 4 of the ndnd states 
with n = 14,16,18 and 20, respectively. We then find all 
two-atom manifolds with and 

within an energy interval of ±5% around Si and denote 
this set of manifolds by Ali. We only retain those man¬ 
ifolds in M-i with orbital angular momentum Ia,Ib < 4. 
For each set of manifolds Aii we choose one atomic sep¬ 
aration Ri such that most quantum numbers within Mi 
give rise to outer turning points Ra and Rb with 


0.8 < 


R^ 

Ra + Rb 


< 1.4. 


( 12 ) 


We find that this regime of weakly overlapping elec¬ 
tron clouds can be adjusted by choosing Ri = 700ao, 
i ?2 = 9 OO 00 , i ?3 = 1200ao and R 4 = 1500ao, and all 
manifolds within Mi that do not obey Eq. (121 are dis¬ 
regarded. After this pre-selection process each set Mi 
typically contains several hundred manifolds, and we ran¬ 
domly select 100 manifolds in Mi that form the set Si. 
Since the total magnetic quantum number is conserved 
by the interaction Hamiltonian V, we confine our analy¬ 
sis to the M = 0 subspace and assign each uaIa mani¬ 
fold in Si a random magnetic quantum number ruA with 
niB = —niA- The distribution of the chosen states with 
respect to the principal quantum numbers ha and ub 
is shown in Fig. [^for all four sets. It follows that each 
set contains a broad distribution of principal quantum 
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FIG. 3. (Color online) Log-log plot of the autoionization 
rates Ea,/ of the randomly chosen states Si as a function of 
SqlQ. Red dots correspond to 5i, green squares show 52, 
blue triangles are for 53 and black diamonds correspond to 
54 . The dashed line at 5q/Q = 10~'* separates regions I and II 
where the autoionization rates behave qualitatively different 
as a function of the overlap. 


numbers where the variation in both ha and ub is larger 
than five. 


V. RESULTS AND DISCUSSION 

The results for the autoionization rate E m of the ran¬ 
domly chosen states in all sets Si are shown in Fig. 
as a function of the overlap Sq/Q, see Eq. Q. There 
are two qualitatively different regions I and H divided by 
the dashed line at Sq/Q = 10“'*. In region I the decay 
rates appear to be independent of the overlap. On the 
contrary, the decay rates increase sharply with Sq/Q in 
region H. In order to understand the physical reason for 
these two regions we perform reference calculations where 
we replace the interaction Hamiltonian V in Eq. 0 by 
its multipole expansion Ume [31] including dipole-dipole, 
dipole-quadrupole and quadrupole-quadrupole interac¬ 
tions. We find that the autoionization rates calculated 
with Ume differ by at most 10% from the values ob¬ 
tained with V for all states with Sq/Q < 10““*. We thus 
conclude that the multipole expansion of the interaction 
Hamiltonian holds if the overlap between the Rydberg 
orbitals is less than 10“^. On the other hand, the results 
obtained by Vme and V differ greatly in region H. While 
all autoionization rates obtained by Ume are smaller than 
3 X 10^s“*, those calculated with V can be several orders 
of magnitude larger for Sq/Q > 10“^. This dramatic 
increase in the autoionization rates with the overlap is 
consistent with the findings in [25]. It can be explained 
physically by noting that the full interaction Hamiltonian 
V allows for direct electron-electron interactions in the 
region where the charge densities overlap, whereas the 
leading term in Ume is the dipole-dipole interaction. 

In the following we analyse the autoionization rates in 
regions I and H in more detail. First, we focus on region 
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FIG. 4. (Color online) (a) Log-linear plot of the autoioniza¬ 
tion rates F m of the randomly chosen states Si in region I as 
a function of energy Em of the initial state. £i is the ce ntra l 
energy of the interval corresponding to set Si (see Sec. IV I, 
and Ry is the Rydberg constant. The dashed line interpolates 
the four mean decay rates (Fm)! [see Eq. @1. (b) Log-linear 
plot of the autoionization rates F m of the randomly chosen 
states Si in region II as a function of SqjQ. The dashed line 
is an exponential ht to the de cay rates Fm in all data sets 
Si with 5q/Q > 10“^ [see Eq. (14l]. In (a) and (b), red dots 
correspond to 5i, green squares show 52, blue triangles are 
for S 3 and black diamonds correspond to 54. 


I and plot all autoionization rates Tm with Sq/Q < 10“^ 
as a function of the energy Em of the initial state |'0 m) as 
shown in Fig.|^a). Within each set Si, the autoionization 
rates show no evident energy dependence. The spread 
in Tm is roughly the same for each set Si and spans 
about two orders of magnitude. However, the lower and 
upper bounds of each set Si depend on energy such that 
the mean decay rates {TM)i become gradually smaller 
by moving from set 5i to 54 , where {TM)i is obtained 
by averaging over all decay rates Tm in 5^ with Sq/Q < 
10“'*. This is illustrated by the dashed line in Fig. 0a) 
interpolating the four mean decay rates {TM)i, 

(Tm) = Ki\EM\/Ryr, (13) 

where k = 2.51xl0*^s“*, 7 = 3.52 and Ry is the Rydberg 
constant. The dominant contribution to the autoioniza¬ 
tion rate in the multipole regime is the dipole-dipole in¬ 
teraction term such that T^ oc dld/.R' where db is 


the dipole matrix element between lips) and I'^b), and dc 
is the dipole matrix element between \iPa) and a Coulomb 
wave. The average principal quantum number n of the 
involved Rydberg states increases from set 5i to set 54 , 
and hence we expect the involved dipole matrix elements 
to increase on average with fi^ m- However, all states in 
a given set Sj a re evaluated at a given atomic separation 
Ri [see Sec. IV with Ri « Ra + Rb oc [see Eq. ([T^]. 
It follows that oc n and hence we expect the full 
autoionization rate to decrease with increasing energy of 
the two-atom state. On the other hand, the large spread 
in Fm within each set Si can be explained with the strong 
dependence of the transition dipole matrix elements on 
the quantum numbers of the initial and bound states. 

Second, we analyze the steep increase of Tm in region 
H. A log-linear plot of the autoionization rates T m in 
region H is shown in Fig. [^b) as a function of Sq/Q. We 
find that Tm increases approximately exponentially for 
Sq/Q > 10“^. This is illustrated by the dashed line in 
Fig.gb) given by 


Tm = rolO° 


(14) 


where the parameters Tq = 3.31 x 10^s“* and a = 15.28 
are obtained by fitting the data points from all sets Si 
with Sq/Q > 10“^ to Eq. (14). The spread of the decay 


rates around the dashed line is roughly three orders of 
magnitude for Sq/Q < 0.15, and reduces to two orders of 
magnitude for Sq/Q > 0.15. In particular, the autoion¬ 
ization rates are apparently independent of the energy of 
the initial state if the overlap exceeds several percent. 

Finally, we note that the overlap of a given state \'4>m) 
is correlated with the symmetry of the energy distribu¬ 
tion between the two atoms. More specifically, we con¬ 
sider the symmetry parameter 


5 = 2- 


.E, 


riAlA 

Em ’ 


(15) 


where En^A is the independent-atom energy of state 
ItpA) and Em is defined in Eq. (10). A value of 5 = 1 
corresponds to a completely symmetric distribution of 
energy Em between atoms A and B, and S decreases 
monotonically with reduced symmetry. Figure shows 
a log-linear plot of S for all sets of states Si as a func¬ 
tion of Sq/Q, demonstrating that symmetry and overlap 
are clearly correlated. This result is relevant for sys¬ 
tems similar to the experimental setup reported in [22j . 
where a gas of cold atoms was excited to ndnd states 
by short laser pulses. This initial state is perfectly sym¬ 
metric with 5=1. However, the interatomic distance of 
some of the atom pairs in the gas will be so small that the 
dipole-dipole interaction couples the initial state to near¬ 
resonant two-atom states with 5 < 1. It follows that 
even if the autoionization rate of the initial ndnd state 
is small for atomic pairs with Sq/Q < 10“^, some of the 
two-atom states involved in the dipole-dipole cascades 
may autoionize much faster because they have 5 < 1 
and hence their overlap can be significantly larger than 
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FIG. 5. (Color online) Log-linear plot of the symmetry S 
[see Eq. @1 of the randomly chosen states Si as a function 
of Sq/Q. Red dots correspond to 5i, green squares show 52, 
blue triangles are for S 3 and black diamonds correspond to 

54. 

for the initial state. A more quantitative analysis of this 
point can be achieved by a simulation of the full quantum 
dynamics starting from an experimentally achievable ini¬ 
tial state and including all coherent couplings between 
two-atom states and their autoionization rates. Such an 
investigation would be an interesting prospect for future 
studies. 


findings hold for other alkali-metal atoms as well. While 
we had to restrict our calculations to relatively small 
principal quantum numbers due to technical reasons, we 
anticipate that qualitatively similar results should hold 
for higher principal quantum numbers as well. Extending 
our current calculations to this regime is subject to fur¬ 
ther investigation. Other possible extensions of our work 
include the calculation of the correct two-electron eigen¬ 
states via full configuration interaction methods [H] , and 
the application of the complex rotation method |50j in 
order to find the energies and widths of the two-electron 
resonances. 

In summary, Rydberg atoms with slightly overlapping 
electron clouds offer fascinating possibilities for future 
theoretical and experimental studies at the boundary be¬ 
tween ultracold atom and molecular physics. In particu¬ 
lar, ultrafast pump-probe laser techniques |I6) allow one 
to resolve processes that are much faster than the au¬ 
toionization rate even if the electron clouds overlap by 
a few percent. In this way autoionization and coherent 
processes in correlated Rydberg electron clouds could be 
measured with unprecedented temporal and spatial res¬ 
olution. Such experiments would represent a paradigm 
shift from mimicking electron-electron interactions with 
ultracold atoms [88ti4^ to actually realizing them. 
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VI. CONCLUSION 


In this paper we present quantum mechanical calcu¬ 
lations for autoionization rates of two nearby Rydberg 
atoms. We consider sets of randomly chosen two-atom 
states and calculate the autoionization rates in lowest or¬ 
der perturbation theory. Since the electron clouds over¬ 
lap only slightly, we neglect exchange corrections to the 
autoionization rate. We find that the autoionization rates 
can be classified via the charge overlap between the two 
states. If the overlap is less than 10“^, the multipole ex¬ 
pansion of the interaction Hamiltonian holds and the au¬ 
toionization rates are relatively small. In particular, they 
decrease on average with increasing energy of the two- 
atom state and can be smaller or comparable to dipole 
transition rates between near-resonant two-atom states. 
It follows that the quantum dynamics in this regime will 
exhibit a rich interplay between coherent transitions and 
autoionization. However, we find that the autoioniza¬ 
tion rates increase dramatically beyond the dipole-dipole 
regime where overlap effects become significant. Our re¬ 
sults show that this regime begins where the overlap ex¬ 
ceeds 10“^, and an approximately exponential increase 
sets in if the overlap is larger than 1%. Our calculations 
were carried out for the specific example of Rubidium 
atoms. However, our classification of the autoionization 
rates in terms of the charge overlap makes no reference 
to the quantum numbers of the initial states or specific 
properties of Rubidium atoms. We thus expect that our 


We thank the National Research Foundation and 
the Ministry of Education of Singapore for sup¬ 
port. The authors would like to acknowledge the 
use of the University of Oxford Advanced Research 
Computing (ARC) facility in carrying out this work 
(http://dx.doi.org/I0.528I/zenodo.22558). The research 
leading to these results has received funding from the 
European Research Council under the European Unions 
Seventh Framework Programme (FP7/2007-2013)/ERC 
Grant Agreement no. 319286 Q-MAC. 


Appendix A: Evaluation of the Coulomb matrix 
element 


Here we outline the evaluation of the matrix element 

M = (Ai) 

entering the autoionization rate in Eq. The operator 
V in Eq. Q is a sum of Coulomb interactions l/|r — r'| 
which whe expand as 


r — r' 


00 I 

l—O m— — l 


47r 


+ 1 rL+i 


(A2) 


where r< = min(r, r') and r> = max(r, U). We truncate 
the sum over angular momenta I in Eq. (|A2[) and omit all 
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terms with Z > 15. Here I corresponds to in the main 
text and determines the amount of angular momentum 
that the Coulomb interaction can transfer between the 
electrons. In order to evaluate the matrix element M 
we expand all wavefunctions in terms of spherical har¬ 
monics [44]. Since we place atom A at the origin, the 
expansion of |' 0 ^) and comprise only a single term, 


Mr) = Rn.iAr)Y/:^e,cj,), (A3a) 

i’Lir) = CEiir)Yr{OA)- (A3b) 


The wavefunctions iPbA) and ipbi'"') Eq. (Al I are cen¬ 
tered at atom B. They are both of the form 


can be found using the orthonormality of Tj ’, 


= J d9d(t)sin6il>i3{r,9,4>)Y^’‘iO,(j)) ■ (A6) 

We represent all radial functions on a grid with up to 
14000 points. The integration region in Eq. (A 6 | is re¬ 


stricted to the solid angle where ippA) takes on non- 
negligible values, and the integral is carried out us¬ 
ing the trapezoidal rule (JS). With the expansions in 
Eqs. (A2|, (A31 and (A5| the evaluation of the matrix 


element M can be reduced to a double integral over the 
radial variables and the remaining integrals reduce to 
Gaunt coefficients [15] . The radial integrals are evaluated 
with the trapezoidal rule [45] , and the Gaunt coefficients 
are defined as 


='lpnplpmp{r - R) (A4) 

with ipnpipmA'^) = and /3 G {b,B}. 

Since R is different from zero we make a general ansatz 
for the expansion of in terms of spherical harmon¬ 

ics, 

^max 

Iq — 0 TTiq — Iq 

where we expressed r in terms of spherical coordinates 
r{r,9A)- We set Lmax = 1000 , and the function Q^^ir) 


GTArr^ = I sin 9Y,A {9, ^Y.^ {0, AAT 

_(T\m 3 / (2^1 -I- 1)(2/2 + 1) 

^ ^ Y 47r(2;3 + 1) 

(Zi, ^2; 0,01^3,0)(Zi, ^2; toi, m2|Z3, —m^). 

(A7) 

The evaluation of Glebsch-Gordan coefficients [35] 
(Zi, ^ 2 ; wi, W 2 IZ 3 , m 3 ) involves the calculation of factori¬ 
als which can be numerically unstable for large values 
of Zi, I 2 and Z 3 if floating point numbers are used. In 
order to circumvent this problem, we generate a library 
of all non-zero Gaunt coefficients with h, I 2 ^ 1000 and 
I 3 < 15 with the software packet MATHEMATIGA [47] . 
The calculation of is implemented in MATLAB [35| . 
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